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ABSTRACT 

Jr  ■  '  - 

-We"  study  the  convergence  of  a  finite  difference  scheme  for  the  Cauchy 
problem  for  the  porous  medium  equation*  u^  =  (u  ,  m  >  1  . 

The  scheme  exhibits  the  following  two  features.  The  first  is  that  it 
employs  a  discretization  of  the  known  interface  condition  for  the  propagation 

7*1  '  J  “  r  '  '  1  *'  4-  »  r  -Z- 

of  the  support  of  the  solution.  We  thus  generate  approximate  interfaces  as 
well  as  an  approximate  solution. 

The  second  feature  is  that  it  contains  a  vanishing  viscosity  term.  This 
term  permits  an  estimate  of  the  form  l(um  1 )  I  <  c/t  . 

_  - .  XX  1  /  R  _  _ 

-  v  We  prove  that  both  the  approximate  solution  and  the  approximate 
interfaces  converge  to  the  correct  ones. 

Finally  error  bounds  for  both  solution  and  free  boundaries  are  proved  in 
terms  of  the  mesh  parameters . 
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SIGNIFICANCE  AND  EXPLANATION 


The  porous  median  equation  represents  a  model  for  the  flow  of  a  gas  in  a 
medium  with  constant  porosity.  The  most  interesting  feature  of  such  a 
physical  phenomenon  is  that  if  initially  the  gas  is  concentrated  in  a  compact 
set,  then  for  later  times  even  though  the  support  expands,  it  remains 
compact.  The  boundary  of  the  region  where  the  gas  is  concentrated  is  the  free 
boundary  of  the  problem.  - - 

In  this  paper  we  construct  numerically  the  solution  of  such  a  problem  and 
also  the  free  boundary. 

We  also  supply  error  bounds  in  terms  of  the  mesh  parameters,  both  for  the 
solution  and  the  interfaces  in  the  sup  norm. 
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AN  INTERFACE  TRACKING  ALGORITHM 
FOR  THE  POROUS  MEDIUM  EQU*mION 

E.  DiBenedetto*  and  D.  C.  Hoff* 


1.  Introduction 


In  this  paper  we  derive  and  analyze  a  finite  difference 
scheme  for  computing  both  the  solution  and  the  interfaces 
for  the  porous  medium  equation  in  one  space  dimension.  We 
demonstrate  that  the  approximate  solutions  and  the  approxi¬ 
mate  interface  curves  converge  to  the  correct  ones,  and 

oo 

obtain  L  bounds  for  the  error  in  terms  of  the  mesh 
parameter. 

Consider  the  laminar  flow  of  a  polytropic  fluid  of 
density  (x,t)  -*  u(x,t)  ,  in  a  porous  medium  which  is 
assumed  to  occupy  the  whole  space,  and  suppose  that  at  time 
t  =  0  the  fluid  is  contained  in  the  slab  £^(0)  <  x  <  tr(0)  . 

The  phenomenon  can  be  modeled  by 

(1.1)  ut  =  (um)xx  ,  (x,t)  £  ST  h  f  x  (0,T] 

0  <  T  <  °°  , 

(1.2)  u( • , 0)  =  uQ(*)  in  R  , 

where  m  >  1  is  a  given  constant,  and  u^  is  a  given  non¬ 
negative  function  such  that  uQ(x)  >0  if  x  £  (C^(0),C  (0)) 
and  u0(x)  =  0  elsewhere.  We  assume  u^  is  continuous  in  1R  . 
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Since  the  problem  is  degenerate,  (1.1)-(1.2)  is  inter 
preted  in  a  weak  sense  and  the  solution  possesses  a  modest 
degree  of  regularity.  Precisely  (x,t)  -►  u(x,t)  is  said 
to  be  a  weak  solution  of  (1.1)-(1.2)  if 

u  €  C(ST)  ;  (um)x  £  L2(St) 

and 

(1.3)  /  u(x,*)<Kx,*)dx|  J  + 

R  l() 


+  f  f  [ - u4>  +  (um)  d>  ldxdt  =  0 

t  R  1  A  A 

ro  K 


for  all  <J>  satisfying 


(1.4) 


<t>  c  H1(St)  Ol°°(St)  and 
x  -*■  ( x , t )  is  compactly  supported  in  R 

uniformly  in  t  , 


and  for  all  intervals  [tQ,t]  -  [ 0  > T ]  • 

The  pressure  v  in  the  fluid  is  connected  to  the 
density  by 


(1.5) 


v  =  u 


m-1 


up  to  some  multiplicative  constant,  and  it  satisfies 


(1.6) 


v.  =  mvv  +  — r  (v  )  ,  in  S„ 

t  xx  m-1  v  x  ’  T 


The  Cauchy  problem  (1.6)-(1.7)  is  also  interpreted 


in  the  weak  sense 


v  €  C(ST)  ;  vxd  L2(St) 


and 

(1.8)  f  v(x,  • ) <J> (x , • )  dx  + 

R  *0 

+  /  /  r “ +  mvvx$x  +  (vx)2<()]dxdt  =  0 

t0  R 

for  all  <j>  satisfying  (1.4),  and  all  intervals  [tQ,t]d 

[0,T]  . 

Existence  and  uniqueness  of  weak  solutions  of  (1.1)- 
(1.2)  was  first  proved  by  Oleinik-Kalashnikov  and  Choui-Lin 
in  [15],  and  the  equivalence  of  (1.3)  and  (1.8)  is  due  to 
Aronson  [2]. 

A  consequence  of  the  degeneracy  is  that  u(«,t)  and 
v(»,t)  are  supported  in  a  finite  interval  [C^(t) ,Cr(t) ] . 

The  curves  (t,C^(t))  ,  (t,cr(t))  ,  which  we  refer  to 

as  the  left  and  right  interfaces  are  Lipschitz  continuous 
and  monotone  decreasing  and  increasing  respectively  (see 
[3]).  The  interface  curves  and  the  pressure  v  are  con¬ 
nected  by  the  Stefan-like  conditions  (see  [3,11]) 


It  should  be  noted  that  conditions  (1*9)  are  not  part 
of  the  original  problem,  but  rather  are  known  to  be  satis¬ 
fied  by  the  unique  solution  of  (1.1)-(1,2). 

Nevertheless  our  algorithm  will  be  based  upon  suitable 
discretization  of  both  (1.6)  and  (1.9). 

We  now  give  a  detailed  description  of  our  algorithm. 

Let  Ax  and  At  denote  increments  in  x  and  t  , 
and  let 


x^  =  kAx  ,  k  £  Z  ;  tn  *  nAt  ,  n  £Mu(0) 


The  approximations  to  v(xk,tn)  ,  e£(tn)  and  Cr (t„) 
will  be  denoted  by  v£  ,  respectively. 

Actually  we  shall  describe  the  computations  only  for 


the  right  hand  interface;  the  computations  near  ^  are 
completely  analogous.  We  therefore  suppress  the  subscript 
and  denote  by  . 

To  start  the  scheme  let  v£  =  vg(x]t)  and  =  c(0)  • 
Next  define 


K(0)  =  max  Ik  : 


x. 


Then  in  analogy  with  (1.9)  we  compute  5 


from  the  equation 


1 


r. 


0 


+ 


m 

m-1 


0 

VK(0) 

s0 


At 


Observe  that  sQ  >  Ax  and  C1  >  ^  . 

Now  given  Cn+1  >  £n  and  v?  for  j  €  Z  ,  we  proceed 
as  follows.  First  define 


(1.10) 


K(n+1) 

K(n+1) 


max  (k  :  xk+1  < 
min  {k  :  xkl  > 


Fig.  1. 

Then  for  K(n+1)  <  k  <  K(n+1)  compute  v£+^  from 


the  finite  difference  equation 
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(1  ...  vk  1'vk  .  ,  n  .  .  vk«r2vk'>vk-l  ,  m  ,vk*rvk-l  .2 

(1.11)  — jj - m(vk  *  O  - —p -  *  srx  ( - jjx -  > 


where  e  >  0  will  be  chosen  later.  Observe  that  we  do  not 
enforce  the  difference  equation  across  the  interface. 

Next  let 


Cl. 12) 


n+1 

s  n+1  "  5  "  xK(n+l) 


and  observe  that 


(1.13) 


Ax  <  S. 


Then  for  xK(n+]j  <  x^  <  £n+*  ,  compute  v£+1  from  the 
linear  interpolation 


(1.14) 


rn+1-x 

_ _ vn+i 

Vi  K(n*1) 


a 

* 

* 


i 


Finally  set  v?  =  0  for  x,  >  en+1  and  compute  sn+2  from 


(1.15) 


_n+2  „n+l  .  in  K( 

r  =  r  +  ■  ■  ~r  - 

m-1  S. 


At 


We  shall  prove  that  v£  >  0  for  all  n  and  k  so  that 
by  (1.15)  £n+2  ^  .  Thus  the  support  of  the  approxi¬ 

mate  solution  increases  monotonically  in  t  . 

In  addition,  the  fact  that  s„  >  Ax  insures  that 

n 

numerical  instabilities  are  avoided  in  the  computations  (1.14) 


and  (1.15). 


Introducing  the  notations 


8  = 


At 


(Ax) 


7  ■  Avk  -  vk.l  -  2vk  +  vk-l  • 


we  can  rewrite  the  difference  scheme  (1.11)  in  the  form 


(1.16) 


n+1  n 


nTi  ii  ,  /•  n  >  ■  ii  ,  mp  (  IV  1  x  iv- x  % 

vk  *  vk  *  m8'vk  +  e)Avk  *  sn:  <■  — 2 —  > 


n  n 

n  .  mB  ,  vk+l"vk-l  ^2 


We  shall  assume  throughout  that 


[Al]  0  <  vQ(x)  <  M  ,  Yx  £  R  ; 

[A2]  |vQ(x)-v0(y) |  <  yq  | x-y |  Vx,y  £  R 


[A3] 


e  is  of  the  order  of  Ax  and 


^  0  m+2  A 

e  >  9  htt:  yoAx  s 


[A4] 


2m8  [M  +  e  +  —j-  yqAx]  <  1  , 


where  M  and  Yq  are  given  positive  constants.  Since 

e  =  0 (Ax)  ,  condition  [A4]  on  8  is  seen  to  be  a  slight 

strengthening  of  the  usual  parabolic  stability  condition. 

We  let  h  denote  the  pair  (Ax, At)  ,  and  we  construct 

approximate  solutions  v*1  and  approximate  interface  curves 
h  h 

and  ?  by  piecewise  linear  interpolation.  Our  results 


may  be  summarized  as  follows 


[I] 

II  vh  -  V II  c 

°°,sT 

<  C (T) (Ax)p 

[II] 

vh  —  V 

X  X 

in  Lq(ST)  , 

for 

all  q  e  [I,00) 

[III] 

,,  ,  h  h, 

^  ^  r 

^£*^r^  lla>,  [0, 

,T]  < 

C (T) (Ax)p/2 

where  p  is  defined  in  terms  of  m  in  Theorem  4.1  below. 
Further  comments  will  be  made  in  Section  6  about  these  rates 
of  convergence,  where  we  present  and  discuss  the  results 
of  some  numerical  experiments. 

The  idea  of  exploiting  an  interface  condition  such  as 
(1.9)  for  computational  purposes  seems  to  have  been  first 
used  by  Huber  [10]  in  connection  with  the  one-phase  Stefan 
problem  (see  also  [7]). 

We  remark  on  the  introduction  of  the  vanishing  viscosity 
e  .  If  e  were  zero,  the  continuous  analog  of  (1.11)- 
(1.15)  would  be  overspecified.  The  artificial  viscosity  e 
thus  seems  to  stabilize  our  finite  difference  scheme.  More 
specifically,  the  presence  of  the  e  allows  us  to  derive  a 
lower  bound  for  v|V  (in  the  sense  of  distributions).  This 

A  A 

in  turn  yields  a  uniform  modulus  of  semicontinuity  for  v” 

h 

and  via  the  interface  condition,  for  5  .  It  is  this  semi¬ 

continuity  which  is  crucial  in  proving  the  convergence  of 
the  approximate  interface  curves,  as  well  as  in  estimating 


the  rate  of  convergence. 

We  briefly  comment  on  related,  at  our  knowledge  available 
results.  In  [9]  Graveleau  and  Jamet  obtained  solutions  of 
the  porous  medium  equation  and  related  equations  by  employ¬ 
ing  a  difference  scheme  similar  to  ours.  However  their 
scheme  is  applied  in  all  of  {t  >  0}  so  that  approximate 
interfaces  are  not  computed.  Moreover  numerical  evidence 
indicates  that  the  supports  of  their  approximate  solutions 
spread  out  too  rapidly  in  time.  Thus  computing  the  inter¬ 
faces  by  "shock  capturing"  seems  to  be  unsatisfactory. 

While  this  paper  was  in  preparation,  Tomoeda  and  Mimura 
f 1 3 ]  informed  us  that  they  have  recently  derived  an  inter¬ 
face  tracking  algorithm  for  the  porous  medium  equation. 
Numerical  evidence  suggests  that  the  approximate  interfaces 
computed  by  their  scheme  are  accurate,  but  they  are  unable 
to  prove  this  result. 

In  addition  their  scheme  is  somewhat  complicated  to 
implement,  since  it  involves  solving  Rieman  problems  for  the 
Burgers  equation  at  each  mesh  point.  Both  their  scheme  and 
ours  suffer  from  the  parabolic  stability  condition 
At  =  0[(Ax)2]  . 

The  paper  is  organized  as  follows.  Section  2  contains 
the  derivation  of  basic  estimates.  Specifically  we  prove 
the  finite  difference  analog  of  the  following  facts,  which 
are  known  to  hold  for  the  exact  solution  v  of  (1.6)-(1.7): 
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(1.17)  0  <  v  <  M 


(1.18) 

(1.19) 

(1.20) 


'V-.s,.  <  Yo  (see  [21) 

I v (x 1 1 2 )  -  v(x,t1)|  <  C|t:2-t1|1/2  (see  [8,12]) 


m-1  1 

vxx  *  m(m+l)  t 


lvxx(‘»t)’vt(*»t)ll,R  <  T  (see  r5]) 


In  Section  3  we  demonstrate  the  convergence  of  the 
approximate  solutions  and  interfaces  to  the  correct  ones 
by  making  use  of  various  compactness  arguments. 

The  error  estimates  are  proved  in  Sections  4  and  5. 

Finally  in  Section  6  we  present  and  discuss  the  results 
of  several  numerical  experiments. 

Throughout  the  paper  we  make  the  convention  that  C 
shall  denote  a  generic  positive  constant  depending  only  on 
m  ,  M  ,  Yq  and  some  specified  time  T  . 


2.  Basic  Estimates 


We  begin  the  analysis  by  establishing  maximum  principles 
for  vj  and  for  the  discrete  space  derivative 


n  _  vk"vk-l 
wk  Ax 


We  assume  throughout  that  the  initial  function  vQ  satisfies 
assumptions  [Al]  and  [A2],  and  that  the  mesh  parameters  e 
and  3  satisfy  [A3]  and  [A4]. 


Lemma  2.1  The  bounds 

(2.1)  0  <  vj  <  m 
and 

(2.2)  Iwjjl  <  Y0 

hold  for  all  k  and  all  n  >  0  . 


Proof.  The  results  hold  for  n  =  0  by  hypothesis. 
Proceeding  by  induction  on  n  ,  we  rewrite  the  last  term 
in  the  difference  equation  (1.16)  as 


12 


& 


>■: 


*, 


- ' 
A 


Ni 


V. 


n  n  n  n 

mAt  ,  wk+l  wk- 1  . ,  vk+l~vk-l  , 
m- 1  ^  2  '  ^  2 Ax 


Rearranging,  (1.16)  thus  becomes 


n+1 


=  [1  -  2mS(vjJ  +  e)]v‘k‘ 


n 


(2.3) 


+  tmB(vk  +  e)  +  i?l 


6Ax  ,  n  +  n*,  n 

^  1  k+1  kJJ  k+1 


r  o  n  .  ->  m  BAx  ,  n  ,  n.  ,  n 

+  [mB(vv  +  e)  -  ^rr  ~r*  (wk+l  wk)]vk-l 


m- . 


Using  the  induction  hypotheses  (2.1)  and  (2.2),  we  have 
that  the  coefficients  of  v£  and  v£+1  in  this  expression 
are  bounded  below  by 


Y0Ax 

1  -  2m6(M  +  e)  and  mB(c  -  ^ (m-1)  ^ 


respectively.  Since  these  quantities  are  nonnegative  by 
[A3]  and  [A4] ,  (2.3)  shows  that  v£+1  is  a  convex  combination 
of  v£  ,  »  and  vk+i  •  This  proves  that  0  <  v£+1  ^  M 

for  k  <  K(n+1)  .  When  k  >  K(n+1)  ,  these  bounds  follow 
from  (1.14). 

We  prove  (2.2)  first  for  k  <  K(n+1)  .  Rewrite  (1.16) 
as  follows: 


(2.4) 


n  A  n 

W  +  w 

n  +  l  n  ^  0 ,  n  A  n  n...  ^  mAt  k  +  1  k-,2 

vk  ■  vk  *  m8(vk  *  c)(Kk*l  -  wk)Ax  *  irr  ( - 3 - ) 


to 
.<  1 


a 


fi  v  ■  f" 
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We  subtract  from  this  the  equation  corresponding  to 
and  divide  by  Ax  .  Using  the  discrete  product  rule 


akbk  "  ak-lbk-l 


ak+ak+l 


<bk  “  bk-i) 


o  . 


we  obtain 


so  that 


bk+bk-l 


(ak  -  ak-i>  • 


n+1 


w. 


-  ,  ,n  r  k  k-1  A  n  ,  n 

-  +  m8  ( - 2 -  +  e)  Aw^ 


wn  -wn 

+  m8Ax  wj  (Jc.4l.?-k:..1) 


wP  +2wn+wn 

.  m  ,k+l  wk  wk-l,  ,.n 

5TT  81  x  ( - * - )twk*l 


-S.i) 


N 


J 

H 

V' 


(2.5) 


w 


n+1 


nx  n 

..n  .  _„,Vvk-l 

wk  *  ™bc — 2 — 


n 


+  e)  Awk 


m  SAx 


m- . 


(wk.i*2",k*wk-i)(wk.i  ■wk-i) 


This  equation  has  the  form 


(2.6) 


where 


wk+1  =  ^  "  2a)wk  +  (a  +  b)wk+i  +  (a  -  b)wk_i  * 


J 
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(2.7) 


a  =  mg  ( — — — * — -  +  c) 


(2.8)  b  ■  Jj  fij£  (wjtl  *  ZmwJ  •  wjl.j) 


By  the  induction  hypotheses  (2.1)  and  (2.2),  a  and  b  satisfy 


(2.9) 


mge  <  a  <  mg(M  +  e) 


(2.10) 


lb  I  < 


i(m+l)  By0Ax 


Thus,  using  the  mesh  conditions  [A3]  and  [A4],  we  obtain 
immediately  that  1  -  2a  and  a  -  lb  I  are  nonnegative. 
Hence  (2.6)  shows  that  is  a  convex  combination  of 

»  and  >  and  so  satisfies  the  bound  (2.2) 

Finally,  we  prove  (2.2)  for  near  the  interface. 

Thus  let  k  =  K(n+1)  and  let  s^  =  cn  -  xk  ,  so  that 


-  vn/s' 
vk/sn 


Using  the  difference  equation  (1.16),  we  then  have  that 
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1  r..n 


sn+l  sn+l 


[vk  +  m6(vk  +  c)(wk+l  '  wk)Ax 


n  n 

.  mAt r  n  ^2  A  ,wk+l  wk^2  ,.n  ^2^ 

*( - 2 - 5  -  (*W  } 


itl  (3-  -  -V  w"  ,  it) 

>_ . ,  *■  n  m-1  k+1  J 


3wn  +wn 

.  ,  n  nv  m  r„A  ,  n  .  At  ,  k+1  kx1 

+  <*W  -  Wk>  •te4*<vk  +  e)  -  ( - 4 - >1 

n+1 


„  _  _ ,  mAt  ,  n 

sn+l  sn  m-1  wk+l 


by  (1.15).  Therefore 


(2.11) 


wn+1  =  wn  +  cfwn 

wk+l  wk+l  clwk 


*!n> 


where 


(2.12) 


c  -  —2—  [BAxfvJ  +  e) 
n+l  K 


7  n  ^  n 
At  ,3<Vl»Wk., 
m-1  ^  ~5  ^  * 


Again,  using  the  induction  hypotheses  (2.1)  and  (2.2),  we 
obtain  that 


(2.13) 


mBAx  V\  _ _  To**, 

- — -  (c  -  — r)  <  c  <  nSCM  +  c  *  — T-) 


the  mesh  conditions  [A3]  and  [A4]  then  imply  that  0  <  c  <  1 


Thus  (2.11)  shows  that  w£*j  is  a  convex  combination  of 
w£+1  and  w£  ,  and  so  satisfies  the  bound  (2.2). 

Finally,  when  k  >  K(n+1)  +  1  ,  w£+*  is  between 
w£*j  and  0  ,  and  so  again  satisfies  (2.2).// 


The  bound  (2.2)  for  w£  ,  together  with  the  interface 
condition  (1.15),  shows  that 


<  C  . 


Since  At  =  0(Ax2)  by  [A4],  it  follows  that  sn  <  2Ax  +  0(Ax2) 
(see  (1.10)  and  (1.12)).  Combining  this  with  (1.13),  we 
therefore  have 

(2.14)  Ax  <  sn  <  3 Ax 

for  small  Ax  .  Actually,  any  upper  bound  on  sn/Ax  will 

suffice  for  our  purposes.  However,  for  the  sake  of  simplicity, 

we  shall  make  use  of  (2.14)  without  mentioning  the  precise 

conditions  on  Ax  which  justify  it. 

In  the  next  lemma  we  establish  a  lower  bound  for  the 

second  spatial  differences  of  v£  .  This  lower  bound  will 

provide  a  uniform  modulus  of  seraicontinui ty  for  w£  and, 

cn+1_cn 

via  (1.15),  for  - — ^  —  as  well.  This  semicontinuity  will 
be  crucial  later  for  obtaining  error  estimates  for  the  approx¬ 
imate  interfaces . 


Lemma  2.2:  Define  K  by 


m-1 

m(m+l)  * 


Then  the  bound 


(2.15) 


.  n  n  ~  n .  n 
Avk  .  vk*i-2vk*vk-i 


Ax 


Ax 


T 


holds  for  all  k  and  all  n  >  0 


_K_ 

t 

n 


Proof:  Denote  the  variable  in  question  by  z£  .  That  is, 


„n  _  Avk  _  wk+l"wk 


Now,  if  0  <  t  <  ,  then 


n 


0 


n,  lwk+ll+lwk!  2y0  K 
k 1  Ax  Ax  t 

n 


We  proceed  by  induction,  assuming  that  (2.15)  holds 
at  time  level  n  ,  and  that 


(2.16) 


KAx 


'n+1 


The  first  case  is  that  in  which  k  <  K(n+1)  -  1  ,  so  that  both 
w£+1  and  w£*J  satisfy  (2.5).  Subtracting  and  dividing  by 
Ax  ,  we  thus  obtain 


=  z“  +  m8( 


„  e)  A2n 


.  mEfVk*1'Vk-l1fAwk*l*Awk1 
*  m6( - Hx - )( - 2 - 5 


*  BAx  r  n  /  ^  t  w  n  .  n»  .  n  i  ✓  _n  „ 

*ST  1 V2  +  (uicn  +  “k>  +  Vil,!w  ‘  z 


Li> 


ffl  |Ax  .  n  ,  0  _n  n  .,11  n  n  n  . 
+  m-1  8  (Zk+l  +  Zk+Zk-l)(wkf2'wk+wk+l-wk-l) 


We  rewrite  the  third  term  on  the  right  as 


mg  Ax  ,  n 


K.l  *  -  zk-l> 


and  the  last  term  on  the  right  as 


FI  T 


The  result  is  that 


(l-2p)Z“  ♦  Cp+q)ZjJ+1  +  (p-q)Z“_1 


(2.17) 


*  r(Z^i  *  ZmZj  *  zJ.1)(Z“n  *  ZZj;  *  Z^.j) 


n  .  „n 


where 


,,n  n 


c  * 
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q  = 


twk+2  +  (4m-1)  ^w*k+i 


n 


n.. 

WJ 


n  i 
"k-i1 


and 


m  At 

r  =  ihT  “5"  * 


We  shall  show  that  z£+1  in  (2.17)  is  an  increasing 
function  of  each  of  the  quantities  z£  ,  ^k-1  *  anc*  ^k+1  * 

Using  Lemma  2.1,  we  have 

,n+l 


-  (l-2p)  +  r[(2m+2)(Za+1+zJ_1)  +  8mZ^] 


n. 


3rj*> 


(l-2p)  +  ^f(2,lH'2)(wH-2  "  Wk-1)  +  (6m“2)  (wk+i“wk)  1 


>  1  -  2m6(M  +  e)  -  ~j  ^  •  16my0 


-  1  -  2m3  [M  +  e  +  ^  YQAx] 


>  0 


by  [A4 ] .  And  similarly, 


fT 

3Z?u 


p  t  q  *  2r  (zjtl  *  (■♦DZj  *  zJ.j] 

PM*  JlK+z  -  “k-l  *  "(wk*l  '  “k)! 


H 


>  mBe  - 


=  mB  [e 


m 


m- 


I 


BAx 


~ir 

3m+l 

(m-l) 


8m>0 

YpAx] 


2m  At 
m- 1  8 Ax 


(2m+2) yq 


>  0 


by  [A3]. 

Thus  is  bounded  below  by  the  right  side  of  (2.17) 

with  Zj  ,  z£_^  ,  and  Z  replaced  by  -K/tn  .  That  is, 

z?*1  >  -  «  .  2fEiL  (JL)2  it 
k  tn  m-1  'V 

K  ri  m(m+l)  K-| 

=  '  11  ‘  m-1  nJ 

_K_  n_-l_  >  _K_  n  K 

'  '  ln  n  '  '  ln  '  '  ^n+1 


as  required. 

There  are  several  cases  to  consider  in  order  to  establish 
the  bound  (2.15)  for  z£+^  when  is  near  .  Now, 

when  z£  +  ^  >  0  ,  (2.15)  is  automatically  satisfied.  We  may 

therefore  assume  that  z£+1  <  0  ;  that  is,  that 


- 2 - 


v 


n+1 

k 


But  this  shows  that  v£+*  is  positive  and  that  v£+1  is  not 
computed  from  the  linear  interpolation  (1.14).  Therefore  it 


c 


.k-\. 

H.  . 

•  N 


.1* 


must  be  that  k  <  K(n  +  1)  .  Since  we  already  dealt  with  the 
case  that  k  <  K(n+1)  -  1  ,  we  may  therefore  assume  that 
k  =  K(n+1)  . 

Thus  and  w£ti  satisfy  (2.6)  and  (2.11)  respec¬ 

tively.  These  equations  may  be  rewritten 


n+1  n  .  _n 

wk+l  "  wk+l  -  cAxZk 


and 


wk+1  =  wk  +  (a+b)Ax  Zl  "  (a-b)Ax  , 


where  a,  b,  and  c  are  as  in  (2.7),  (2.8),  and  (2.12) 
Subtracting  and  dividing  by  Ax  ,  we  thus  obtain 


i 


(2.18)  Z*j+1  -  (l-a-b-c)zjj  ♦  (a-b)zj  j 


We  checked  in  the  proof  of  Lemma  2.1  that  a- lb  I  >  0  . 
And,  using  (2.9),  (2.10),  and  (2.13),  we  have  that 


1-a-b-c  >  1  -  mB(M  +  e)  -  - mB(M  +  e  +  ~ _) 

m—  I  2  m— 1 


m+3 


=  1  -  2mB[M  +  e  +  4 YqAx] 


which  is  nonnegative  by  [A4] 
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Finally,  using  (2.9),  (2.10),  (2.13),  and  (2.14),  we 


have  that  the  sum  of  the  coefficients  in  (2.18)  is 


v  Ax 

l-2b-c  <  1  *  BVx  -  f  (e  -  ^r)  . 


Using  the  definition  of  K  ,  this  bound  may  be  rewritten  as 


1  .  6Y0AX  fl  1  ,  mBe 
1  +  K  [1  +  3li+l) ]  T" 


And  using  condition  [A3] ,  we  then  find  that 


l-2b-c  <  1  - 


20Yq^x 


_  2Y0  At 

'  1  'Tax 


On  the  other  hand,  we  have  from  (2.16)  that  tn+1  >  KAx/2yq  , 


so  that 


1  <  2y0  At 

imr  ^  nr  37  • 


Therefore 


1 - 2b-c  <  1  - 


and  (2.18)  shows  that 


1  _  n 

n+I  "  iuT 


_K_  n 
*n  n57 
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Remark;  The  constant  K  =  (m-l)/m(m+l)  is  the  best  possible 
as  it  can  be  verified  on  the  Barenblatt-Pattle  solution  [6,16] 
(see  also  [5]). 


We  can  improve  the  bound  (2.15)  by  imposing  additional 
regularity  conditions  on  vn  . 


Corollary  2.3;  a)  If  vQ  is  a  concave  function,  then 


(2.19) 


for  all  k  and  all  n  >  0  . 

b)  If  there  is  a  constant  CQ  such  that 


(2.20) 


vn(x+h) -2vn (x) +vn (x-h) 


>  -  C, 


holds  for  all  x  and  all  h  >  0  ,  then 


(2.21) 


^  -  Cr 


for  all  k  and  all  n  >  0  . 


Proof:  We  showed  in  the  proof  of  Lemma  2.2  that,  in 
all  cases,  Av£+Vax^  is  an  increasing  function  of  Av^/Ax^ 
and  Av^^/Ax^  .  The  bounds  (2.19)  and  (2.21)  then  follow 
easily  from  (2.17)  and  (2.18).  // 


* 
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Next,  we  obtain  a  bound  for  the  discrete  time  derivative 
near  the  interface. 


Lemma  2.4:  There  is  a  constant  C  such  that 


n+1  n 

!vx  -vk 


holds  for  n  >  0  and  k  >  K(n+1)  . 


Proof:  We  use  the  symbol  0  (•)  to  denote  dependence  on 

mesh  parameters.  Let  k  >  K  =  K(m+1)  and  let  s'  =  cn  -  x,,  . 

n  k 

Then  from  (1.14)  we  have  that 

vn+1-vn  cn+1-x  tn-x 

vk  vk  _  1  xk  n+1  5  xk  n. 


n+1  n  ,n+l 

VK  ”VK  1  ^  ~ Xk 

°cd  i-V-  +  ^  <-r— 

n+1 


t  ”*ki  ni 
~V~)  VK 

n 


The  first  term  on  the  right  can  be  estimated  by  using  the 
difference  equation,  (1.16): 

vr1-vK  „  WK+ 1-WK 

-V1  ■  *  O'1’  ’  “f1’ 


since  v^+e  =  O(Ax)  .  And  the  second  term  on  the  right  is 


O(Ax) 


[^n+0(At) -xk] s^- [ s'+0(At) 1 (Cn-xk) 


n  bn+l 


=  0(1)  .  // 
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Lemmas  2.2  and  2.4,  together  with  the  difference  equa¬ 
tion  (1.11),  now  imply  the  following  bound  for  the  discrete 
time  derivative. 

Corollary  2.5:  There  is  a  constant  C  such  that 
n+1  n 

v,  -v,  . 

(2.22)  ■— -  C(1  +  -L) 

Ln 

holds  for  all  k  and  all  n  >  0  . 


In  the  next  lemma  we  use  the  one-sided  bounds  (2.15) 

and  (2.22)  to  derive  L  estimates  for  Av?/AxZ  and 
..n+1  n 


-v. 


At 


Lemma  2.6:  a)  For  a  given  T  >  0  there  is  a  constant  C 
such  that 


Av 


n 


Ax 


Ax 


Ax  ^  C(1 


holds  for  tn  <  T  . 

b)  If  the  initial  function  Vq  satisfies  the  hypothesis 
(2.20)  ,  then 


Av 


n 


Ax 


Ax 


.n+1 


v. 


n 


At 


Ax  <  C 
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c)  And  if  Vq  both  satisfies  (2.20)  and  is  concave,  then 


Avk 

„n+l  n 
vk  -vk 

Ax7 

9 

At 

Proof:  From  Lemma  2.2  we  have  that 


2K 

t 

n 


We  multiply  by  Ax  and  sum  over  k  .  Since  v£  is  zero 
outside  an  interval  of  length  C(1  +  t  )  ,  we  obtain  that 


Ax  <  (1  +  t  ) 

ln 


The  other  bounds  in  a)  and  b)  are  proved  similarly,  c) 
follows  from  Corollary  2.3,  the  difference  equation  (1.11), 
and  Lemma  2.4.  // 

In  the  final  lemma  of  this  section  we  establish  the 
Holder  continuity  in  time  of  the  sequence  (v£)  • 


Lemma  2.7:  Let  T  >  0  be  given.  Then  there  is  a  constant 
C  such  that  the  inequality 

!VU-Vkl  <  cdV1.'1'2  *  s  (c  *  e'^IV'i/72 

holds  for  t  and  t  in  [0,T]  and  for  all  k  . 
n  m  ’ 


-  V"  1  f  ■  ■ 


1 


a 


I 


V, 

i 


f: 


p 


B 


t: 
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Proof:  Fix  a  point  (x,  ,t  )  and  let  t  >  t  be 

kq  nQ  nj  nQ 

given.  Let  Q  be  the  rectangle 


Q  *  [x,  -  p,x,  +  p]  x  [t  ,t  ] 

R.Q  Kq  IIq  11^ 


where  p  is  a  multiple  of  Ax  to  be  chosen  later.  And 
define  the  quantities 


H  *  max  |vP  -v.°|  , 

ng<n<n^  K0  K0 


(2.23)  c  *  2m(M  +  e)  +  yQp  , 


and 


uk  -  <  ■  vk°  -  V  -  7  t(xk-xk„32  ♦  1 


We  shall  show  that  <  0  for  Cx]t»tn)  £  Q  by  induction 
on  n  .  When  n  *  nQ  and  |x-XjJ  ^  P  »  we  have,  using 
Lemma  2.1,  that 


_  Hrt  Ha 

Uk  vk  ~  vk  ~  Y0P  K  0 


For  the  induction  step  we  consider  the  following 
three  cases : 

Jxk“xk  I  *  p  »  k  >  K(n+1)  ,  and  |xv-xv  I  <  P  with  k*%K(m+l) 


k  Ak, 


If 
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In  the  first  of  these,  we  have  that 


u?*1  <  tvj-x*1)  *  cvr1-^ 

K  K  KQ  Kq  Kq 


YqP  -  H  , 


which  is  nonpositive  by  (2.2)  and  the  definition  of  H  . 

In  the  second  case,  we  have  v£+*  ^  Y0sn+1  ^  3  YgAx  *  so  t^iat 


Uk+1  <  3  Y0Ax  "  Y0P  » 


which  is  nonpositive  provided  that 


(2.24) 


p  >  3  Ax 


For  the  third  case,  we  employ  the  linearized  difference 
operator  L  ,  defined,  for  a  given  sequence  z£  ,  by 

zn+l_zn  AZn 

,  7n+l  Ak  „r„n  A  ^  AZ,k 

Lzk - At - *Cvk  *  e)  77 

Vn  -vn  Zn  -Zn 

m  fVk+l  vk-l./k+l  Vl, 
HTT[ - Jftx - n - TEx - }  * 


Applying  L  to  u£  and  using  (1.11),  we  find  that 


-T  t-c  +  2m(vk+e)  +  STT( 


n  n 

m  ,VVl  ,Vl‘Vl 


xk  ^ 
0 


so  that,  by  (2.1)  and  (2.2), 


Ill 
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LUk+1  ("c  1  2m(M+e)  + 


m-. 


Y0PJ 


=  0 


by  the  definition  of  c  ,  (2.23). 

On  the  other  hand,  we  can  rewrite  the  inequality 
LU”+1  <0  in  the  form 

Ujj+1  <  [1  -  2mB(vJ+e)]u£ 


n 


+  [mB(vk+e) 


'k+l"vk-l 


)]U 


k+1 


n , 


+  [m6(vk+e)  + 


B 

J 


m 

m- 


(v 


n 


k+1  vk-l'Juk-l 


The  coefficients  of  U?  on  the  right  hand  side  of  this 


inequality  are  exactly  the  same  as  those  of  va  in  equation 

(2.3).  And  we  showed  in  the  proof  of  Lemma  2.1  that  these 

,  1 

coefficients  are  nonnegative.  We  therefore  have  that 


is  a  convex  combination  of  U 


k-1  *  uk  »  and  Uk+1  » 


U 


and 


so  is  nonpositive  by  the  induction  hypothesis. 


Setting  k  =  k. 


in  the  result  u£  <  0  ,  we  thus  obtain 


that 


n, 


Y0P 


*“!■ 

p 


J 
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where  s  *  t  -  t  .  In  a  similar  way,  we  can  establish 
n,  nn 

*  ®  nQ  rj 

the  same  inequality  for  v,  -  v,  .  Taking  the  maximum  over 

K0  0 

n  6  [nQ,^]  ,  we  thus  obtain 


(2.25)  H  <  y0P  +  ^  s 

P 

We  shall  choose  p  so  that  <  1/2  .  Specifically, 
p  should  satisfy 


p^  +  3Ax  <  p  <  p^  +  4Ax 


where  p^  is  the  larger  root  of  the  quadratic  equation 


p2  -  2cs  *  p2  -  YqSp  -  4m(M+e)s  =  0 

1/2  1/2 

An  easy  computation  shows  that  p^  =  0(s  +  s  )  =  0(s  ) 

for  t  <  T  .  And  p  >  3Ax  ,  as  required  by  (2.24).  Since 
ni 

<  1/2  ,  (2.25)  becomes 

PZ 

II  <  2yqP  <  2Y0(p1  +  4Ax)  *  C(s1/2  +  Ax) 


In  particular, 


n  i  n_ 

lvk  “vk  I  *  cdtn  -tn 

k0  K0  nl  n( 


1/2 


+  Ax) 


// 


We  remark  that  the  proof  of  the  above  lemma  is  the  discrete 
version  of  an  argument  given  by  Kruzkov,  [12]  and  Gilding  [8], 
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3.  Convergence  of  the  Approximate  Solutions 


Let  h  denote  a  pair  (Ax, At)  whose  elements  satisfy 
the  mesh  conditions  [A3]  and  [A4 ] .  We  define  approximate 
interface  curves  t  i— * ►  ^(t)  and  t  *— *-  c^(t)  by  piecewise 
linear  interpolation:  for  t  <  t  <  tn+i  * 


cjct) 


'n  +  K(n)  (t-t  ) 


where  cn  ,  v£  ,  and  S  are  as  in  Sections  1  and  2;  and 
similarly  for  cj(t)  .  The  estimate  (2.2)  for  vJJ(n)/sn  then 
shows  that  the  nets  U^(t)}  and  u£(t)}  are  uniformly 
Lipschitz  and  uniformly  bounded  in  finite  time. 

t  We  construct  approxi- 

tn+l..  mate  solutions  v^(x,t)  in 

an  analogous  way,  as  follows. 

Tn 

t  .  \S _ 111  I  f  T?  and  S?  are  the 

n  k  k 

triangles  in  Fig.  3.1,  then 


Fig.  3.1 


vh(x,t)  =  v£  +  (x-xk)wjj  +  1  ♦  »  U»t)  €  t£  ; 


vh(x,t)  =  vk+i  +  (x_xk+l)wk+l  +  (t-tn+l)ok  *  ^  Sk 


It  follows  immediately  from  Lemmas  2.1,  2.2,  2.6,  and  2.7 


that  the  functions  v  satisfy 


(3.1)  0  <  vh(x,t) 

(3.2) 

»£  <Y 
^5T  T0 

a  2  h 

(3.3) 

— V(x,t) 

^ _ “ 

and 


a.e. 


3x 


^r(x.t) 


for  t  >  0  are 


finite  measures  in  R  with  mass 


I 

R 


.2  h 

3  V 


3x 


■l  fr"  < 


(3.4)  |vh(x,t+s)  -  vh(x,t)|  <  C(T)  s1/2  ,  0  <  t  <  t+s  <  T 


Throughout  this  section  we  fix  a  time  T  and  a  rectangle 
Q  =  [a,bj  x  [0,T]  ,  where  [a,b]  is  large  enough  to  contain 
the  supports  of  v^(*,t)  for  all  h  and  all  t  £  [0,T]  . 

(See  (1.15)  and  (2.2)  ) 

h  li  u 

The  properties  of  ^  ,  and  vn  described  above 

insure  that,  for  every  sequence  h.  tending  to  0  subject 


to  the  mesh  conditions  [A3]  and  [A4],  there  is  a  subsequence, 
which  we  index  simply  by  h  ,  for  which 

v*1  — *•  v*  uniformly  in  Q  , 

?! uniformly  in  [0,T]  , 

and  v£  — ►  v*  in  L2(Q)  ,  weakly  . 

Our  goal  in  this  section  will  be  to  prove  that  v*, 

5*  ,  and  coincide  with  the  exact  solution  and  inter¬ 

face  curves  for  the  problem  (1.8).  Actually,  the  convergence 
li  li  h 

of  v  ,  e"  ,  and  also  follows  from  the  error  bounds 

which  we  shall  derive  later  in  Sections  4  and  5.  However, 
the  arguments  of  the  present  section  are  much  more  direct. 
Moreover,  we  obtain  here  the  convergence  of  v^  in  L^(Q) 
for  all  p  <  00  .  As  a  byproduct  of  these  arguments,  we 
thus  obtain  in  addition  a  constructive  proof  of  the  exis¬ 
tence  and  regularity  properties  of  the  solution  of  (1.8). 

We  begin  by  showing  that  v^  -*■  w*  strongly  in  L^(Q) 
and  that,  in  fact,  w*  =  v*  . 

Lemma  3.1:  For  any  t  >  0  the  net  (v^(*,t)}  is  pre- 

1  A 

compact  in  L*[a,b]  . 


Proof:  The  proof  consists  of  estimating  the  L1 

difference  between  v^(*,t}  and  spatial  translates. 

Given  t  >  0  ,  choose  n  so  that  t  <  t  <  t  . 

7  n  n+ 1 

Then  when  h  is  sufficiently  small,  0  <  t-At  <  t  . 

If  we  take  p  =  £Ax  where  l  is  a  positive  integer, 
then  it  is  easy  to  see  that 


(3.5) 


/  |  vjj(x+p,t)  -  v£(x,t)  |  dx 


<  c  I  wk+^"wk  I  Ax  +  Zklwk+£"wk  lAx-) 


The  first  of  these  sums  is  bounded  by 


s  £W  'Xti'  ix2 

EkEj=k+l  Ax 


i  n  n 
w.  -w. 


i.  U.  ■ 


j  k=j -l I  Ax 


li  AX2 


|Av?| 


-  p  Z.  - i-  Ax  < 


by  Lemma  2.2.  Dealing  with  the  second  sum  in  (3.5)  in  a 
similar  way,  we  find  that 


f  |v£(x+p,t)-v£(x,t)  |dx  < 
R 


(3.6) 
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holds  when  p/Ax  is  a  positive  integer. 

When  p^  <  Ax  ,  | v^Cx+p^ ,t) -v^(x,t) |  will  be  zero 

except  when  x  is  within  p^  of  the  nonhorizontal  sides 
of  the  triangles  t£  and  s£  .  Thus 


f  | vjj(x+p1,t)-vjj(x,t)  |dx 


<  Cp1(£k  |wj+1-wj|  +  £k  lwk+J~wk  X|) 


Pi 

But  this  expression  is  ^  times  the  right  side  of  (3.5) 
with  l  taken  to  be  1  .  The  computations  we  made  above 
therefore  show  that 


l_  1.  C  P  -m 

f  |vJ(x*Pl,t)-vJ(x,t)|dx  <  ^ 

Combining  this  with  (3.6),  we  see  that  (3.6)  now  holds  for 

all  p  >  0  .  The  conclusion  of  the  lemma  now  follows  from  [1],  // 


Lemma  3.1  thus  shows  that  {v£>  has  strong  L*  limit 
points.  In  the  next  lemma,  we  prove  that  these  limit  points 
can  be  identified  as  the  derivatives  of  limit  points  of  {v*1}  . 

L 

Lemma  3.2:  Let  vn  denote  a  sequence  of  approximate  solu¬ 
tions  which  converge  to  v*  uniformly  in  Q  .  Then 
h  n 

v"  -*■  v*  in  LP(Q)  for  every  p  £.  [  1 , «>)  . 

^  A 


ll 

Proof:  By  Lemma  3.1,  every  subsequence  of  {v£(*,t)} 
has  a  subsequence  which  converges  in  L^(R)  .  Thus  let 
{h'l  <=  { h }  and  let  v^  (*,t)  converge  to  a  function 
£(x)  in  L*(R)  .  We  shall  show  that  £(x)  =  v*(*,t)  a.e. 
First,  if  \p  6  H^fajh)  ,  then 

f  [vjj  (x,t)-ip  (x,t)]  [vh*  (x,t)-i|;(x,t)  ]  =  0 
a  x  x 

We  take 

<J>  =  n<J>  +  (l-n)v*(*  ,t) 

1  h* 

in  this  relation,  where  4>  £H0(a,b)  .  Since  v  -*•  v* 

h»  1 

uniformly  and  vx  (*,t)  -*■  £  in  L  we  obtain,  by  letting 
h  -*■  0  ,  that 

b 

n  /  [C-n4>  -  (l-n)v*(»  ,t)  ]  [<j>-v*(*  ,t)  ]dx  =  0 
a  x  x 

Dividing  by  n  and  letting  n  ■*  0  ,  we  then  find  that 
b 

f  U-v*(*  ,t)  ]  [4>-v*0  ,t)  ]dx  =  0 
a  x 

for  all  <f>  £  H^(a,b)  .  This  shows  that  E,  =  v*(*,t)  a.e. 
Thus  every  subsequence  of  (vn(*,t)}  has  in  turn  a  sub¬ 
sequence  which  converges  to  v£0,t)  in  L^CR)  .  And 
therefore  the  entire  sequence  converges  to  v*(*,t)  in 


for  every  h 


9 


we 


L^OR)  .  Finally,  since  llv^ll^  ^  Yo 
have  that 


ivi-v;' 


P.Q 


T 

C(p)[/ 

0 


lvx(-,t)-vj(.,t) 


1,R 


dt] 


1/p 


-  0 


by  the  dominated  convergence  theorem.  // 

We  remark  that  the  proof  of  Lemma  3.2  is  an  adaptation 
of  an  argument  given  by  Minty  in  [14], 

The  next  theorem  contains  the  main  results  of  this  section 


Theorem  3.3:  Let  v  ,  ^  ,  and  cr  denote  the  exact  solu¬ 

tion  and  interface  curves  for  the  problem  (1.8).  Then 


(3.7) 


v  uniformly  in  Q 


(3.8)  v£  — *■  vx  in  LP(Q)  ,  p  <  00  , 

and 

(3.9)  — *■  Cr  uniformly  in  [0,T]  • 


We  prove  (3.7)  by  showing  that  limits  of  converging 
sequences  from  {v^1}  satisfy  the  weak  equation  (1.8),  and 

so  agree  with  its  solution  v  ,  which  is  known  to  be  unique. 

H  h 

Thus  let  {v  }  denote  such  a  sequence  and  let  vn  ■*  v* 
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9 


uniformly  in  Q  ,  so  that  v^O,l)  -*■  v*(*,t)  in  L1(R) 
for  every  t  >  0  .  It  will  be  sufficient  to  show  that 


N 


Fi 


/  v*(x,  *)4»(x,  •)  |T  dx  + 
R  *1 


(3.10) 


2 

/  /  [ - v* 0  +mv*v*<f>  +  ^ (v*)24>]dxdt  =  0 

ri  L  X  X  III*  X  Jv 


T1 R 


for  all  C  functions  4>  satisfying  (1.4)  and  for 


0  <  T.  <  T,  <  T  . 


Given  such  a  function  <J>  ,  let  <J>£  =  ‘Hx]c»tn)  and 
consider  the  quantity 


(3.11) 


„  .  n+1  n  .  n  n  n 

N_-l  v.  -v,  .  Av,  v.  ..-v,  ,  0  „ 

„  2  r^rk  k  ,  n,  k  m,  k+1  k-ls/T  ,n. 

En-N,  {Ek'— —  -  '  1^1  <  Mx  >  ’  VX,“  ' 

1  Ax 


for  appropriate  and  N2  .  Now,  the  expression  in  brackets 

vanishes  for  K(n+1)  <  k  <  K(n+1)  .  And  for  other 

vrx 

values  of  k  ,  we  have  that  - ^ -  =0(1)  by  Lemma  2.4, 


.  n  n  n 

n  Avk  wk+l"wk 

and  (v£+e) — 7  -  0(Ax) - t -  =  0(1)  by  Lemma  2.1.  Thus  the 

K  Ax^ 


Ax 


quantity  (3.11)  approaches  0  as  h  ■>  0  . 

On  the  other  hand,  we  can  sum  by  parts  in  (3.11)  and 
match  the  resulting  terms  with  the  corresponding  integrals 
in  (3.10).  We  shall  carry  out  the  details  only  for  the  most 
complicated  term.  Using  Lemma  2.2,  we  may  rewrite  the  second 
two  terms  in  (3.11)  as  follows: 


39 


n  n  n  n 

„„  r,  n,  x  ,n  Wk+1  Wk  ,  1  ,Wk+l  kx2.ii!.  . ^ 
-m  EE  [<vk+eMk  —  +  ^  ( - 5  >  ^JAxAt 


n  (vk+£)<-(vnk-l+E)Cl  1  n  2  n, 

m  EE  [vP  {— ! - -  - - —}  -  Z7T  (*£)♦?  AxAt  +  O(Ax) 

k  Ax 


m-1  v  k'  TkJ 


(3.12) 


n  n 

V  —V 

rr.  r  n  r  k  k-1  x  n  /  _n  , \ 

m  11  K  {  Ax" *k  +  (V’-  1+0 


|,%n 

pk  vk-l 


k-1 


4x  >  *  ST  <wk>2+fcItoAe  +  0(4x> 


EE  (w^)2^AxAt  +  m  EE  v^wj<t>x(x”, tR) AxAt  +  O(Ax) 


for  some  x£  £.  [xk _  ]  .  We  shall  show  that  the  second 

sum  here  converges  to  the  second  term  on  the  right  in  (3.10). 
First  note  that,  since  v^(x,t)  =  w£  on  T£uSk  1  , 

wP AxAt  *  ff  v*Vx,t)dxdt 

K  __  _  i  X 

i i cn “ 1 

Tkubk 

And,  since  <j>  is  smooth  and  is  Lipschitz  in  x  and 

Holder  continuous  in  t  , 


= 

ff  V^V^4>  dxdt  +  0(Ax+At^2)  AxAt 

Tn  U  Sn" ^ 
k  U  bk 

Therefore  the  second  sum  in  (3.12)  is 
T2 

/  /  mvhv^4>  dxdt  +  0(Ax)  , 

Tj  R  x  x 


which  approaches 
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/  /  mv*v*4>  dxdt 

Tj  1R  X 

as  h  -+  0  ,  since  v*1  -*■  v*  in  l”  and  ->  v*  in  L1  , 

X  X 

L  1* 

with  vn  and  uniformly  bounded.  The  other  terms  in 

(3.11)  are  handled  in  a  similar  manner.  Thus  v*  satisfies 
(3.10)  for  all  appropriate  test  functions,  and  so  coincides 
with  the  unique  solution  v  of  the  problem  (1.8). 

The  proof  of  (3.7)  is  based  upon  the  following  techni¬ 
cal  lemma,  which  will  be  used  again  in  Section  4  for  the 
derivation  of  error  bounds  for  the  approximate  interface 
curves . 

Lemma  3.4:  Let  { }  be  a  subsequence  such  that  -*■  (,* 
uniformly  in  [0,T]  as  h  ->  0  .  Then  for  every  t  >  0  and 
for  any  positive  numbers  6  and  n  , 

t+n  -»  p  ,2 

(3.13)  /  vUJ(s)  -  6,s)ds  >  {  *  U*(t+n)  -  C*Ct)]  -  }  , 

dc* 

for  0<t,  t  +  n^T.  Moreover,  if  exists  and  is 

positive  at  t  ,  then  there  are  positive  numbers  and  C 

such  that 

(3.14)  vUJ(t)  -  6,t)  >  C6 
holds  for  0  <  <S  <  <sfi  . 
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Proof:  Let  p  ,  q  ,  and  N  denote  the  largest  integers 

in  6/Ax  ,  n/At  ,  and  t/At  respectively.  Then 


(3.15) 


N+q-1 

E 

n-N 


VK(n)-p 


At 


N+q-1 

E 

n=N 


K(n) 

E  w  Ax]At  . 

j=K(n)-p+l  J 


using  Lemma  2.2  and  the  definition  (1.15)  of  ,  we  have 
that 


K(n)  K(n)  K(n-l)  Av* 

-  E  w  Ax  ■  E  [-w  ,  \  +  — j  Ax]  Ax 

j“K(n)-p+l  ^  j=K(n)-p+l  n  Ax 


K(n)  ^(t  .J-tNt  >  !  nAx 

Z  l  -r—  — - - r-  n - ±-  22L  iAx 

m  _  1  At  m+1  t  J 

j=K(n)-p+l  n 


m-1 


Slk  [^r(tn+l)~^r^tn)  _  JL  (pAx)! 

1  At  F  mtl  t 


m 


■]  Ax 


Substituting  this  into  (3.15)  and  discarding  the  nonnegative 


term  v 


n 


K(n) 


we  obtain 


(3.16) 


N+q-1 


n=N 


K(n)-p 


N 


Now  (3.13)  follows  by  letting  h  -*■  0  and  using  the  uniform 

h  h 

convergence  of  vn  and  s"  . 
dc* 

If  (t)  >  0  ,  then  there  is  a  positive  number  p 


such  that 
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cj(t+n)  >  s£(t)  +  pn 


h 

Proof  of  3.9:  Let  denote  any  subsequence  con¬ 

verging  to  a  curve  t*  uniformly  in  [0,T].  We  shall  show 
that  C*(t)  =  Cr(t)  for  every  T  .  First  observe  that, 
since  v^(x,t)  =  0  for  x  >  t^(t)  ,  v(x,t)  must  be  0 
for  x  >  C*(t)  thus  Cp(t)  <  C*(t)  . 

Now  suppose  that  5r  <  c*  on  (t,t+n)  with  ^r(t)  =  z. 
Then  since  cr  is  increasing,  there  must  be  a  time 
t  €  (t,t+n)  at  which  d^*/dt  exists  and  is  positive.  But 
then  (3.14)  shows  that  v(s*(t)  -  6,t)  is  positive  for  small 
6  .  However,  this  implies  that  Sr(t)  >  t*(t)  ,  which  is 
a  contradiction.  Therefore  there  is  no  maximal  time  t  for 
which  C*(t)  =  ^r(t)  for  0  <  t  <  t  .  And  since  sr  and 
agree  at  t  *  0  ,  they  agree  for  all  t  .  Similar  argu¬ 
ments  hold  for  ^(t)  •  // 


i-i  * 
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4.  Error  Bounds  for  the  Approximate  Solution  and  Interface  Curves 


In  this  section  we  prove  the  following  theorem. 


Theorem  4.1:  Fix  T  >  0  .  Then  there  is  a  constant  C  such 


that,  for  0  <  t  <  T  , 


(4.1) 


l(*,t)-v(.,t)llo%R  <  C  min  [  (A.X“J  )P*  ,t1/2+Ax] 


(4.2)  |  Ch  (t)  -  C  (t)  |  <  C  ti7T(AxCt|logAx|)ZCp+3) 


where  c  is  either  or  e  .  Here 

c  r 


,  1  <  m  <  2 


,  2  <  m 


,  1  <  m  <  2 


m+1  ,  2  <  m 


We  remark  that,  if  the  initial  data  Vq  satisfies  the 
hypothesis  (2.20),  then  the  term  |logAx|  may  be  omitted 
from  the  bounds  in  (4.1)  and  (4.2). 
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lu 


The  proof  of  Theorem  4.1  will  be  given  in  a  sequence 
of  lemmas.  First  we  introduce  the  weak  truncation  error 
associated  with  an  approximate  solution.  If  <J>  is  a  smooth 
function  satisfying  (1.4),  define 


(4.3) 


2 

J  ( v  t  <j) ,  t  x  ,  1 2 )  =  /  /[vt4>  ♦  mvvx'5)x  +  v^4>]dxdt 


Thus  v  is  a  weak  solution  of  (1.6)  if  and  only  if 
J(v,<J>,tj  ,t2)  =  0  for  all  4>  and  all  intervals  (*1**2^  . 

The  weak  truncation  error  associated  with  the  approximate 

h  h 

solution  v  is  then  the  functional  J (vn+e ,  • , t-^ , t2)  .  We 

have  the  following  estimate  for  J  . 


Lemma  4.2:  Let  f  satisfy  (1.4)  and  assume  that  f^  and 

f  are  in  L°°(IRx[0,T]  with  f(x,T)  =  0  for  all  x  .  Then 

x  2-m 

*  1  f  h  m-T  . 

if  <J>  =  (v  +e)  f  , 


(4.4) 


where 


| J(vh+e,4>,6,T)  |  <C  lllflll  Axa|log6 


fill  =  Ilf  II  +  Ilf,  II  +  Ilf  II 

00  £  co  X  00 


The  proof  of  Lemma  4.2  is  quite  technical  and  lengthy, 
We  therefore  postpone  it  to  Section  5. 
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Next,  we  define  functions  u  and  u  by  the  relations 


u  *  v 


h  ,  li  ,m-. 
u  =  (v  +  e) 


In  the  following  lemma  we  exploit  the  above  estimate  for 
J  to  obtain  a  bound  for  u^1  -  u  in  Lm+*  . 


Lemma  4.3:  There  is  a  constant  C  such  that 


X  b  X 

/  /  |uh-u|m+1dxdt  +  /r/CCuh)m  -  Um)v  dt ] 2dx  < 
0  a  1R  0 

(4.5) 

m 

<  C  [L(u,uh)  Axa  |  log  Ax  |  +  (Ax)5577] 


where 


L(u,uh)  -  +  l/((uh)“-um)xdtI.jRx[0_I] 


Proof:  Let  <p  and  f  be  as  in  Lemma  4.2.  Then  after 
integrating  by  parts  in  (4.3)  and  substituting,  we  obtain 


J (vh+e , <}) ,  <5  ,T)  =  (m-1)  f  f  [ujf  +  (uh)™f  ]dxdt 

6  R  1  x  x 


(4.6) 


*  (m-1)  [-/  u  (x,6)f(x,6)dx 


-/  /  (uhft-(uh)^fx)dxdt] 
5  R  '' 
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We  shall  replace  6  in  these  integrals  by  0  .  The 
resulting  error  in  the  second  integral  on  the  right  will  be 
no  more  than  CHI  £  III  6  .  And  to  estimate  the  change  in  the 
first  integral,  we  use  the  fact  that  uh  is  Holder  con¬ 
tinuous  in  t  ,  which  may  be  established  as  follows.  When 
m  £  (1,2]  we  have  from  Lemma  2.7  that 

|uh(x,t2)-uh(x,t1) |  <  C  |vh(x,t2)-vh(x,t1)| 

<  c  |t2-t1|1/2 

And  when  m  >  2  , 

|uh(x,t2)-uh(x,t1) T'1  <  iuh(x,t2)m"1-uh(x,t1)m' 

=  |vh(x,t2)-vh(x,t1) |  ^c  |t2-t1|1/2 

Thus  u*1  is  Holder  continuous  in  t  with  exponent  a/2  , 
and 


|  /[uh(x,6)f(x,6)-uh(x,0)f(x,0)]dx|  <  C  lllflll  6  a/2 
]R 


We  thus  obtain  from  (4.6)  and  Lemma  4.2  that 


/  uh(x,0)f(x,0)dx  +  f  f  [uhf  -(uh)mf  ] dxdt 
1R  0  R  z  x  x 


Clll  fill  (Axa  |  log6  |  +  6a/2) 


Now  subtract  from  this  the  weak  form  of  the  original  equation 


for  u  ,  (1.3),  and  take  6  =  Ax^  .  The  result  is  that 


/  I  [(uh-u) f  -  ((uh)m-um)  fldxdt 
OF  r  x  X 


(4.7) 


<  C  III  f  III  Axa  !  logAx 


We  have  used  here  the  fact  that 


!lu(*,0)-uh(*,0)Ha>^  <  CAx“ 


which  follows  directly  from  the  definitions  of  u  and  u*1 
and  from  the  Lipschitz  continuity  of  vQ  . 

In  (4.7)  choose 


f(x,t)  =  / 
T 


[(uh)m- 


um] (x,s)ds  + 


(T-t) 


It  is  clear  that  f  satisfies  (1.4)  and  the  conditions  of 
Lemma  4.2.  Observe  also  that 


T  h 

/  /  (u  ~ u)  f , 
0  F 


T  111 

>  »  f  h  v  r  f  h\  m  m  m- 
/  /  (u  -u) [ (u  )  -u  -e 

0  F 


] 


T  b 
>  f  f 
0  a 


uh-u| 


m+1 


-  C  e 


am 


and  that 
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r  r  h\  m  in  i  r  1  9  r  ^  h  m  m  >«  /  \  »  t  2 

[(u  )  -u  ]xfx  =  J  ^  ((u  )  "u  )x(x,s)ds] 


so  that 


/  /  [(uh)m-um]  £  =  -i  /[/((uh)m-um)  (x,s)ds]2dx 

OR  x  R  0 


Now  (4.5)  follows  by  making  these  substitutions  into  (4.7).  // 


The  next  lemma  contains  the  corresponding  estimate  for 
h 

the  error  in  vn  . 


Lemma  4.4:  There  is  a  constant  C  such  that 


/  f  |vh-v)Pdxdt  <  C[L(u,uh) Axa | logAx]  +  Axam] 
0  a 


where  L  is  as  in  Lemma  4.3  and  a  and  p  are  as  in 
Theorem  4.1. 


Proof:  When  m  >  2  ,  we  have 


(vn+e)-v|  = 


Uh  \  m  ■  1  m  *  1 1  i  « .m- 1  ■  h  ■ 

u  )  -u  j  <  (m-l)M  [u  - u | 


And  if  1  <  m  <  2  ,  then 


1  1 


(vh+e)-v|"^T  <  |  (vh+e)"rT-vm 


uh-ul 


Since  e  =  O(Ax)  ,  we  have  in  either  case  that 


|vh-v|p  <  C(|uh-u|m+1  +  Axp) 

The  conclusion  then  follows  from  Lemma  4.3.  // 

00  h 

In  order  to  deduce  the  L  bound  for  v-v  from  the 
above  Lp  bound,  we  shall  require  the  following  interpola¬ 
tion  inequality. 


Lemma  4,5:  Let  z  be  a  continuous  function  on  ST  =  IRx[0,T] 
with  support  z(*,t)  <z-  [a,b]  for  t  £.  [0,T]  .  Assume  also 
that 


a)  zx  £  L°°(St)  , 

and 

b)  zt(»,t)  €  L1^)  with 


lztC*,t)ll1 


Let  p  6  [1,°°)  and  Z  >  0  be  given.  Then  there  is  a  con 
stant  C  independent  of  Z  such  that  the  inequality 

C«)P^  T 

holds  for  0  <  t  <  T-£  .  Here 


F(z)  = 


[  sup  llz(*,s)ll1'  *  +  Cn  II  z 


>,ST  t<s<T 


On 

X 


Proof:  For  t  >  0  we  have 


z(x,t)|p+1  <  (p+1)  /  |z(x,t)|p|z  (x,t)|dx 


(4.8) 


C  "  Zx"“,0  >^P 


where 


y(t)  =  • 


Differentiating  the  definition  of  y(t)  ,  we  obtain 


(4.9)  y(t)P_1 1^2.|  «  CU(.,t)|P-l  lzt(-,t)«ljR 

Combining  (4.8)  and  (4.9),  we  thus  obtain 

P(P-1)  P-1 

/(t)13"1!^!  <  Cy(t)  P+1  lzx|P^  II  z t  ( - ,  t; 


so  that 


El*  E^i 

.P^T.dy i  <  c  ii,  .P7? 


(4.10)  y(t)p  <  c  lz  |P  i  lz.(  - ,t) I, 
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Now  let  xft)  be  the  cut-off  function 

0  <  t  <  T-i 

T-t  <  t  <  T 

Then  for  q  >  1  and  0  <  t  <  T -l  ,  we  have 
y(t)q  =  (y  x) q  I  t 

<  q  /  (xy)q"1( lxt \y  *  xlg^l) 

/  (yq  -  yq'1l^l) 

■I 

Now  choose 


and  use  (4.10)  to  estimate  the  second  term  in  the  above 
integrand.  The  result  is  that 

p(p+3)  T  J±.  2zk 

y(t)  11+1  <£[/y(t)Pdt]  sup  [  yCs)1*1  +  Hz  llP+* 

'  0  t<s<T  X  ’  r 


llzt(*»s)lll^ 


52 


Using  hypothesis  b) ,  this  may  be  rewritten  as 


(4.11) 


p(p+3)  2p  P-1 

y(0  ptl  <h"2"p,s  1  +  c„izji^  ] 


T  t<s<T 


pJR  Ox  »,s 


Now  by  (4.8)  we  have  that 


(4.12) 


!z(* ,t) 


1  _E_ 

<  cilzxiP^  yU)^ 


The  result  then  follows  by  substituting  (4.11)  into  (4.12).  // 


We  can  now  prove  the  first  half  of  Theorem  4.1  by  applying 
the  interpolation  inequality  of  Lemma  4.5  to  the  error 

bound  in  Lemma  4.4  as  follows. 


Proof  of  (4.1):  It  is  clear  that  v  -v  satisfies  the 
first  hypothesis  of  Lemma  4.5.  In  addition,  hypothesis  b) 
is  satisfied  because  of  (3.3)  and  the  results  of  [5]. 
(Alternatively,  we  can  observe  that  the  conclusion  of  Lemma 
4.5  remains  valid  for  limits  of  functions  which  satisfy 
hypotheses  a)  and  b)  uniformly.  And  v  is  such  a  function 
by  the  results  of  Section  3.)  Replacing  T-f  and  T  by 
T  and  T+l  respectively,  we  therefore  obtain  that 


Iv  (*,t)-v(*,t) 


IV  -VI 


where  S.j.  =  Rx[o,T+l]  and  p  £  (1,“)  .  But  if  we  choose  p 
as  in  the  statement  of  Theorem  4.1,  then  Lemma  4.4  shows  that 


I 

Therefore 


II vh - v II P  c,  <  C  Axa |  logAx | 
P  j 


llvh(*,t)-v(-,t)lloo  R 


On  the  other  hand,  the  Holder  continuity  of  vh  and  v  in 
time  shows  that 


llvh(*,t)-v(*,t)ll^]R  <  C(t1/2  +  Ax) 

for  any  t  .  The  estimate  (4.1)  follows  from  these  last 
two  inequalities.  // 


Next,  we  deduce  the  bound  (4.2)  for  the  error  in  the 
approximate  interfaces  from  the  above  bound  (4.1)  for  the 
error  in  v^  .  Again  we  drop  the  subscripts  and  denote 
sjj  and  Cr  by  £h  and  £  . 

Proof  of  (4.2):  First  we  refine  the  result  in  Lemma  3.4 


Divide  (3.13)  by  n  and  let  n  -*■  0  .  The  result  is  that 
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(4.13) 


vU(t)-y,t)  [Jc(t)  -  Ti^TF]  a.e. 


A  similar  bound  holds  for  vn  :  taking  q  =  1  in  (3.16), 


we  have  that 


vK(n)-j  38  TT  !(jAx>  '  $rh~}  ■ 


h 


Since  v  is  Lipschitz  in  x  and  Holder  continuous  in  t  , 
and  since  5*1  is  piecewise  linear,  we  can  conclude  that 


(4.14) 


vhUh(t)-y,t)  >  [^hCt)y  -  (Xl)t]  "  CAx  a*e< 


Now  let  {v1,v2>  =  {v,vn}  and  let  ^  and  c2  be  the 
corresponding  interface  curves.  Assume  that,  at  time  t  , 

0  <  ?2(t)-c1(t)  =  y  .  Then  we  have,  using  either  (4.13)  or 

(4.14),  that 


lvh(*  ,t)-v(*  ,t)  l^R  >  l(v2-v1)(C1(t),t) 


(4.15) 


=  v2(c1(t),t)  =  v2 (s2 (t) -y ,t) 


m-1  •  • 

>  — ^2'^  (m+1)  t  1  '  CAx 


Now  let  F(t)  =  ]  and  let  E(t)  be  the 

bound  in  (4.1)  for  II  v^(  •  ,t) -v(  •  ,t)  II  .  (4.15)  then  shows 

that 

-  T5STTt  F(t)  4  E(t)  - 

where  we  have  subsumed  the  CAx  term  into  E(t)  .  Integrating 
we  thus  obtain  that,  for  any  6  >  0  , 

2  .2  t  _  2 

(4.16)  F(t)  <  t55571  [6  55171  F ( 6 )  +  /  s  55571  E  (s)  ds  ] 

6 

Now,  since 


F ( 6 )  -  Uh(6)-5(<S)]2 

<  ( | Ch(6)-Ch(0) |  +  | C (0) - c ( 6)  |)2 

<  C62 


the  first  term  on  the  right  of  (4.16)  approaches  0  as  6-^0 
In  addition,  the  integrand  in  the  second  term  of  (4.16)  is 
bounded  by 


1 

pT? 


(Ax®  |  logAx | ) 
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A  short  computation  shows  that  the  exponent  of  s  here  is 
greater  than  -1  ,  so  that  this  term  is  integrable  on  [0,t]  . 

We  may  therefore  conclude  that 

1  1 

Uh(t)-s(t)|  =  F(t) 1/2  <  Ct^  (Axa|logAx|)2(p+^  .  // 


Estimate  for  the  weak  truncation  error 


This  section  is  devoted  to  the  proof  of  Lemma  4.2.  Since 
the  arguments  are  technically  involved,  even  though  simple, 
we  begin  by  listing  a  few  facts  to  which  we  will  refer  system¬ 
atically.  In  the  estimates  to  follow  we  will  use  repeatedly 
the  bounds  established  in  Section  2  without  specific  mention. 

Let  and  s£  be  the  triangles  in  Fig.  3.1,  and  let 

v*1  be  the  piecewise  linear  interpolated  of  v£  introduced 
in  Section  3.  In  ,  v*1  can  be  written  in  any  one  of  the 

following  equivalent  forms 


v"  *  wk*i(x-xk)  *  0kti(t-V 


(S.1) 


vh(x,t)  =  |  v£+1  +  wk+].(x-xk+1)  +  °k+i(t‘tn) 


■  vk+l  +  wk+l(x'xk+l)  +  °k+l(t_tn+l) 


,n  . 


and  in  Sk  it  can  be  written  in  any  one  of  the  equivalent  forms 


vk  +  wk!^x-xk5  +  ak(t-V 


(5.2) 


vh(x't}  =  +  wktl(x'V  +  °k^t-tn+l^ 


v?*i  +  wfthx-x,^)  +  a.n(t-t  A.) 
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Let  $  e  H^CS^,)  be  a  test 

L 

denote  with  <j>  the  piecewise 


,n 


values  <J>V  =  4>(xwtn)  and  set 


,n  An 

,n  _  V*k-1 
’''k  ZZ - 


function  satisfying  (1.4), 
linear  interpolation  of  the 


$ 


n 

k 


It 


in  rl 


<P  can  be  written  in  any  of  the  equivalent  forms 


(S.3) 


4>h(x,t) 


f  *"  *  *Lltx-xk)  * 

■  *Li  *  *kn(x-xkn)  *  *i*iCt'V 
.  *k!j  *  *k.i(x-xk.1)  ♦  ♦kn^'Vi5 


and  in 


(5.4) 


<J>h(x,t) 


^k  +  ^k+1Cx-xk)  +  *?(*-*„) 

*k+1  +  ^k+1Cx"xk)  +  $k^t-tn+1) 

Ci  +  ^k+1tx-xk+l)  +  *k(t-tn+l) 


Remark  From  (5.1)  —  (5.2)  it  follows  that  v^1  can  be  written 
as  the  value  of  v*1  at  any  one  of  the  corners  of  (s£ 

respectively) ,  plus  terms  of  the  order  of  Ax  .  An  analogous 
fact  holds  for  <t>^  . 


5.1  The  basic  identity 


Consider  the  quantity 


(5.5) 


J(vh  +  e,4>h,6,T)  = 


. ,  ,  hAh  t  ,  h  .  h,h  ,  m(m-2)  ,  h->2,h,, 
ff  ( vt4>  +  m(v  +e)vx<(>x  ♦  mi ^  --  (vx)  4>  Idxdt 


where  6  is  a  fixed  positive  number.  We  assume  for  simplicity 
that  T  ■  (N+l) At  and  6  *  nQAt  for  two  positive  integers 
nQ  <  N  ,  and  calculate  the  various  parts  of  (5.5)  as  follows. 


ff  v^<f>hdxdt  * 

«  t 


E  Z  {  ff  vVdxdt  +  ff  vVdxdt) 

nft  Z  Tn  cn 

0  lk  bk 


l  Z  *  Z  Z 

n0Z  nQ  Z 


o^l(Ax)2At 


J2  *  m  ff  { (vh+e) vx<fix  +  (vx)24>h}  dxdt  * 
ST,  6 


m  ff  v^  [  (vh+e)  4»h]  dxdt 

n  A  a 

&T,<5 


Z  Z  {  ff  v£+1[(vh+e)<}>h]xdxdT  +  ff  [(vh+e)  d>b]xdxdx  } 
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On  t£  we  have 


[  (vh+e)  4>h] . 


,  s  ,n  ^  n  .n 

fVE)Vi  *  * 


Therefore 


*  vAi|(x->iki  *  lx-xk.i)]  * 


*  *  "wil'f'-'n'  • 


f!  w",1[(vh*E)*h]xdxdT  =  *  wJ^jJAxit  ♦ 

*  *  wk*i*k.i!w"»l  ^(dt)2  . 


By  direct  calculation 


[Av] 


wk+l[Cvk+e)^k+l  +  wk+l^k+l^  =  '(vk+e)  2  ^k  + 

(Ax) 


wn  -wn 
wk+l  wk  .n 


*n  rt>n 

+  f„n+ci  "k+1  "k  .n  .  ,11..  n  ^k+l'^k  A  ,  n  .2. 

^vk  ^  £x  ^k+1  ^vk  wk  Ax  ^wk+l^  ^k+1 


n 


-  ■ ( vV+e) 


[Av] 


n 


(Ax) 


k  ,n  .  1  r  /*,  n  .  .  n  .n  ,  nx  ,  n,n, 

7  +  k  *  lx  ((vkU  +  E)wk*l»W  -  Cvk+e)wk4>k> 


This  implies  that 


Z  I  //  w“vl[Cv“+e)r]xdxdt 


II  Q  Z  (pH 
k 


1  N  n  [Av]"  n 

■ill  (V, +e)  - ft-  4>,  AxAt  + 


n0  Z 


*1  Z  I  +  w£$J]Ax(At)' 


n0  Z 


Analogous  calculations  on  give 


I  r.  ff  w{J^[(vh+e)4>h]xdxdt 

no  z  s£ 


1  ?  r  (vn*l  .  [Avlk  1  n+1 

7  „£„  z  (Vk  1 1SF  k 


1  V  r  „n+l  r„n,n+l  +  wn+1*ni  AvrAt- v 

F  }  l  wk+l  t°k^k+l  wk+l$klAx(At) 


n0  Z 


Substituting  these  calculations  in  the  expression  of  J 
we  obtain 


l  l  m(vj?+e) 


'k  xn 


v1  z 


f  Ax') ' 
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m  nn  [Av]k  n0  N+l  [Av]k  N+l 

-  T  At  Z  Ax{  (v,  °  +  e)  - °  +  (v£+1+e)  - K-  4»”+1} 

1  Z  k  (Ax)Z  K  K  (AX)Z 


,  id  r»  r*  n  /  n  -  n  ix .  n  •»  «  / . .  \  2 

*  F  „E  i  Wk  *  Vk)4x(At) 
n0  2 


m  r  -  n+lf  n.n+1  .  n+l.n^A 

”  (T  Z  Z  wk+lCcVk+l  wk+1$k)Ax(At) 

n0  2 

We  finally  transform  the  remaining  integral  in  J(v^+e ,<|>\ 6 ,T) 


Ciii) 


J,  =  -  -2tr  //  (vh)2d>h  dxdt 
j  ra-  x  q  a 

1,6 


-r  i  n  //  c^+1)24>hdxdt  +  //  c^tj)2<i>hdxdt> 

nn  Z  r**n  cn 

u  rk  bk 


Using  (5. 3) -(5. 4),  by  standard  calculations  we  obtain 


J3  * 


i?T  *  I  K)2*"Ax4t  ♦ 


n0  +  1  2 


*  ttctt  e  l  («J)2<CAx)24t  - 


STTFTT  L  z  UwJ*1)2  -  (wJ)2]^Ax(At)2 

no  2 


At  i  <|tcw“#>2*“°  *  (w^v^1]  - 


lrr  n(k  2  . n0  .  ,  N+1n2  .N+l-i  i a 
"  ?■[  (wi,  )  4*1,  +  (w.  )  <^v.  ]  }Ax  • 


Fuwk  J  "'k 


We  rewrite  the  first  summand  in  as  follows 


n  A  n 
w,_ .  -,+W. 


n  A  n 
w,  .  ,+w,. 


n-*2xn  _  ,  k+1  k^2^n  A  r,  n.2  f  k+1  ks2,  n 
»k)  4>k  ■  ( 5 J  4>v  +  LCwv)  -  ( 7 J  J<J>V 


and  combine  the  expressions  of  J.  ,  i  =  1,2,3  so  obtained, 
as  parts  of  J(vh+e,4)\6,T)  to  deduce  the  following  basic 
identity 


(5.6) 


J(vh*e,*h,6,T)  - 


-I  E  { 

V1  Z 


vr-kn 


,  ».  ,  [Av]k  -  , 

(Ax) 


n  n 

m  ,Vk+l~Vk-l.  2,  ,  n 


L>  X 


+  At  Z  ak04)kAx  - 

z 


mA,  nn  l  Ay]  n  ju.1  [Av],  ^ 

2—  l  {(v\)  - 4>.°  +  (v^+e)  - O^JAx  + 

2  Z  k  (Ax)2  k  k  (Ax) i  k 


*  ^  Z  Z  [oM*}  -  a%J]AxAt  - 
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mAt  _  v  n+lr  n,n+l  .  n+l.nlA 

-5-  z  i  Vi["k*w  *  "kti<fk)4xAt  * 


no  z 


IN  j 

1  l  (w£)  i^AxAt  + 
n0  2 


mAt  n  r-.  n  r  n  .  n  ,  n « n  -1 A  » 

*  -5-  s  I  Wk  *  VkJAx4t 

n0  2 


mAt  ~  r  r  n+l>.2 

nJ  *  [(wk  1 
n0 


(w£) 2]$"AxAt 


mAt  v  r/  no.2.no^  ,  N+lN2,N+ln  lr,  \2,n0  ,  ,  MtL2,N+llu  . 

‘  2^  *  f(wk  >  ♦k  +<Wk  >  ^k  >  *k  +(wk  >  *k  i}Ax  + 


*  Trrrr  A,  l  rrr?  [3w£  *  »k.i]<AxAt  • 


^rr  „;,i  i  7^ 


10 

=  l  H. 

i  =  l  1 


From  now  on  we  will  select  test  functions  <j>  of  the  form 


<{>  =  (vh+e)HfrT  f 


m  >  1 


where  f  satisfies  (1.4).  We  will  estimate  the  , 

i  =  1,2,.  ..,10  in  (5.6)  in  terms  of  f  ,  f  ,  f  .  For 

A  C 


notational  simplicity  set 


x  *  A  t  °°,  S, 


Ilf  II  ■  Ifl  *  If  I  *  If  I 


We  start  by  making  elementary  estimates  of  <f>  and  its  first 
derivatives  in  terms  of  f  . 

Denoting  by  <j>^  the  piecewise  linear  interpolated  of 
$  and  recalling  that  e  is  of  the  order  of  Ax  we  have 


(5.7) 


ll<|>hll  < 

T  oo 


C  Ilf II 


if  1  <  m  <  2 


C  II  f  II  CAx)m" A  if  m  >  2 


Also 


, n  _ .n  2 - m  2  - m 

,, n  1  r  r  n  x  m- 1  rn  ,  n  .  m- 1  _ 

*k  ■  — E - Ix{<Ve)  fk  -  (vk-l*e)  f 


2-m  2-m  2-m 

1  rr/.n.^ra-l  ,  n  ,  •,  ,n  .  ,  n  .  xm-l  k  k-1 

Sx^V0  -  (v^+e)  ]}fk+(\.l+£)  — ~ 


<  I ^ttI  C e ) m” 


I  w”  I  Ilf  II  ♦  (vC_1+e)HrTllf 


(b) 


,nfl  .n 
♦l.  -K 


2-m 


_  k  k  ^  2-m i,  ft.  vm-l 

— —  4  trl<VE> 


-  1 


2-m 

I ^ k  I II  f  H  +  (v^+e)m-1llftll  , 


where 


n 


v.  = 


f  n  n  , 
min  {vk  ;  vk-l} 


n 


v,_  = 


min  {vj+1  ,  v?} 


From  (a) -(b)  we  deduce 


(5.8) 


l*kl 


c  (v^+e)  “ 1 II  £11  +  C  II  f 


if  1  <  m  <  2 


2-m 
,  m- 


C  [  (v^+e)  -1 II  f  II  +  Hfvll](Ax),,,"J-  ,  if  m  >  2 
k  x 


(5.9)  \*l\  <  < 


C(v5J+e)_1|oJ|  II  f  II  +  C  It ft  I 


if  1  <  m  <  2 


2-m 
,  m- . 


[  C  [  (v^J+e) -1 1  0]J  [  II  f  II  +  II  f  t  II  ]  (Ax)J,,_x  if  m  >  2 


5.2  An  auxiliary  lemma 


Lemma  5.1:  For  all  test  functions  <p  of  the  form 


2-m 

.  /  h  \  m  —  T  - 

<P  =  (v  +e)  f  , 


the  following  estimate  holds 


| J(vh+e,4>  ,S,T)  |  <  C  III  £111  (Ax)  a  |  log  6 


9 


where 

a  =  min  (1  ; 

Proof:  We  estimate  the  on  the  right  side  o£  (5.6) 

separately.  At  points  where  the  difference  equation  (1.11) 
holds,  the  summand  in  vanishes,  and  by  virtue  of  Lemma 

2.4,  the  sum  extended  over  the  remaining  (n,k)  is  of  the 
order  of  Ax  . 

Therefore 

|HX  I  <  C  (Ax)»<t>Sgo  <  C  lllflll  (Ax)a 
We  have  also  easily 

|H2|  +  |H3|  <  C  lllflll  (Ax)“ 

In  estimating  we  only  consider  the  term  whose 

summand  is  •  The  estimate  for  the  term  whose  sum¬ 
mand  is  is  analogous,  and  in  fact  simpler. 

We  have 

N  N 

Ax  Z  Z  ojJipjJ^AxAt  =  Ax  Z  Z  ak^k+  lAxAt 

n0  Z  n0  (k<R(n+l) ] u  [k>K(n+l) ] 


ra 


+  C(~«  )  H  f  11  (Ax) 

v  +  e 


Now  observe  that  =  v1*  +  Ax  =  v1]^,*  +  a ”  At  +  w 

K  k+1  k+1  k+1  k+1 


n  n  n  n+1 

vv-v^  -v _ 

L  K  k+1  j  „n  k+1  k+1 

where  w  =  - j- -  ,  and  a  ^ - 77 - 

k+1  Ax  k+1 


.  By  the 


estimates  in  Section  2 


Iv’L  <  c  Ax 

k+1  k+1 


and  therefore 

vk+£  Ax 

(  -"I'  )  <  1  +  C  ^  <  C  . 

v _ +e 

k+1 

(21 

Consequently  for  1  we  have 

(21  a  N  [Av]k 

|<UJ|  <  C  lllflll  (Ax)  E  E  - ^  AxAt 

d  nQ  Z  (Ax) ^ 

By  Lemma  2.6 

N  [Av]^  N  , 

E  E  - £  AxAt  <  C  E  (1  +  i-)At 

nQ  Z  (Ax)  nQ  n 

<  C(l  +  I  log  6  |  ) 


so  that 
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|  <  C  lllflll  (Ax)a  |  log  6  | 


estimate  Ik^  j!  by  performing  a  discrete  integration  by  parts 

„  n  ,  n  ,n+l  ,rri-l 

N  w,  .  1+w.  .  .--(fr, 

-S-  Ax  Z  Z  (Jctlk) 2  (  k+l-k— >  AxAt  ^ 

m-I  t  fix 

n0  K(n+1) <k<K(n+l) 


N  n  +  1 

<  C  Ax  Z  l  - y  j<f>,  I  AxAt  +  C  Axil II  < 

nQ  Z  (Ax) 


<  C  lllflll  (Ax) a  |  log  6 1  . 


Combining  these  estimates  we  obtain 


|  H4  |  <  C  lllflll  (Ax) a  |  log  A  |  . 


By  similar  calculations  involving  the  use  of  Lemma  2.6  we  deduce 


| H5  |  +  | Ily  |  <  C  lllflll  (Ax)a  |  log  6 1 


In  estimating  we  first  integrate  by  parts  (discrete 

integration)  and  use  the  techniques  above  to  obtain 

|  H6  |  <  C  lllflll  (Ax) a  |  log  .5 1  . 
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Analogous  techniques  give  the  desired  estimate  for  Hg  ,  Hg  , 
.  The  proof  is  complete. 

Corollary  5.2.  If  -L  <  vqXx  ^  0  *  ^or  some  Positive 
constant  L  ,  then 


1  J(vh+e,<j>h,6  ,T)  |  <  C  lllflll  (Ax)a  ;  V«  >  0 


5.3.  Proof  of  Lemma  4.2 


Since  J(vh+e  ,T)  -  J(vh+e,<j>-<J>h,6,T)  +  J  (vh+e ,  4>h,  <S  ,T) 
and  the  last  term  has  been  estimated  in  Lemma  5.1,  we  have 
only  to  estimate  J(vh+e ,<J>-<j>h,6 ,T)  .  To  this  end  we  will  need 
the  following  preliminary  fact 


Lemma  5.3.  There  exists  a  constant  C  independent  of  n  , 
k  ,  Ax  ,  4>  such  that 


(Ax'At)"1  ff  Cd>h-4>) dxdt  <  C  lllflll  (Ax) 01  [l+(vj+e)_1] 


min  v 

rpn 


where 
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i 


n 


Remark  An  analogous  statement  holds  for 


Proof  of  Lemma  5.3 


I 


If  C 4>  -<J>)dxdt  = 


nn 


IS  ( 0 dxdt  + 
Tn  K 

K 


i 


+  si  |>J(x-xk)  +  <Dk(t-tn)  ]dxdt  =  I 

Tl 


For  t  using  (5.8)-(5.9)  we  have 


|I(2)|  <  C  lllflll  (Ax) a  [l+(vJJ+e)"1]AxAt 


~n, -1 


Estimating  I 


(1) 


« 


-  //  t  (4>^-0 Cx, tn)  ]  [ <J> ( x , t )  -<f>(x,t 

„,n  ^ 


{  II  <f>  II  Ax  +  II4).II  At}  AxAt 

A  mil  L  mil 

"•Tk  ">Tk 


and  the  desired  estimate  follows  from  (5.8)-(5.9). 


* 


Corollary  5.4:  For  all  n  ,  k  ,  if  vk  =  max  v*1 


(v£+e)  , 

A'x'At  ff  (♦  -*)dxdt  <  C  lllflll  (Ax) 
Tn 
‘k 


r" 

a 


(i)  +  j (2) 


n)  ] } dxdt  < 
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Remark : 


//  (<j>h-<)>) dxdt  <  C  IHflll  (Ax)a  . 

Tn 

Me 

Corollary  5.4  holds  if  T?  is  replaced  by  S1} 


h  hi 

We  are  now  in  the  position  to  estimate  J(v  +  e,4>-4>  ,6,T) 


(5.10) 


J (vh+e  ,<j>-<f>h ,  6 , T)  =  ff  (v^C4>-4>h)  +  m(vh+e)  (vh)x(<i>-<J>h)x  + 

ST ,  6 


+  Cm(vx^  -  ~TCv^2)(<l»"<l>h)>dxdt  = 


=  l  E  [oj 
nQ  Z  1 


m  wf  ,+wjJ  -  , 

^(—4 - )  ]  "  (*-<r)dxdt  + 


+  l  l  [a 

n0  Z 


J"  ff  (<j>-<})h)  dxdt  + 

jn~  x  _ 


i  E  l  ^3wk+l+wk+2^wk+l-wk  +  2^  ^  )dxdt  + 

n0  Z  Tf 


+  j z  1  (3wkti+wk+1)(wk+1’wkli)  n  (^-4>h)dxdt  + 

no  z  s£ 


+  m  //  [(vh+e)  ( vh )  C 4> - 4>h )  +  (vx)  2  (<f-4>h)  ] dxdt  =  I  P- 

i , 5  1-1 


Estimate  of  | |  +  |  P 


P1  = 


N 
E 

IT 


n  .  n 

Wi,.  ,+Wi 


e  [ok  "  ^rr(  — ^>2]  ff  (♦-♦h)dxdt  + 

0  K(n+l)<k<ac(n+l)  S? 


N 

+  E 


n0  [k<K(tttl)  ]  U  [k>K(n+l)  ] 


[o“  "  ^l(Wkfr~^)21  ff  (<M>Vxdt 


=  P*  +  PJ* 


For  P**  by  the  estimates  in  Section  2  and  Corollary  5.4 
we  obtain  easily 

|  PJ*  j  <  C  lllflll  (Ax) a 

In  estimating  P*  we  use  the  difference  equation  (1.11). 

ff  (4>-4>h)  dxdt  |  } 
sk 

<  C  infill  (Ax)°  I  log  6  I  . 

The  estimate  of  | P?  |  is  analogous. 


P1 


N 

C  I 
n„ 


lAv]J 


(Ax)' 


AxAt  { 


(vk+e) 

AxAt 


Estimate  of  | P 3 1  +  j  | : 


|P3|  +  |P4|  <  C E 
n0 


E 

Z 


[Avlk 

(Ax)2 


AxAt{ 


Ax 

AxAt 


ff  (4>— 4>h)  dxdt  |  }  . 

s£ut£ 


By  Corollary  5.4  and  Lemma  2.6 

|P3I  ♦  |P4I  <  C  Ilf  III  (Ax) a  |  log  6  |  . 


Estimate  of  P,.  : 


N 


Pr  « 


m  E  Ziff  [  (vh+e)  ♦  (vV  (<|>-4>h)  ]  dxdt  ♦ 

__  -»  _  A  A  A 

nn  Z  Tn 

u  Ak 


+  ff  [(vh+e)v*C*-<!>h)  +  (v^)2(<J>-<|>h)]dxdt} 

Sn 


=  oM  ♦ 


+  Q 


We  estimate  Q^1-) 


in  x  over  each  T?  . 


by  performing  an  integration  by  parts 


(cn,,tk> 


<t„-xkn> 


We  have 


If  (vh+e)v^(<Ji-(j)h)xdxdt 
n 
Tk 


tn+l  Xk+1 

I  I 

t  .Ax  . 

”  xk+4t*t_tn) 


(vh+e)  (v^)(^-$tl)^dxdt 


“  1  u  u  v, 

/  (vh+e)v^C0-<l>n)  (t,xk+1)dt  - 


/  (vh+e) v^(<{>-<J>h)  (t,Xj,  +  ||  - 


ff  (vV(<|>-<|>h)dxdt 

Tk 


Therefore 


Q(1)  =  m  I  l  {  n/+1Cvh+e)w?+1(0-0h)(t,xk+1)dt  - 

n0  2  tn 

tn+l 

-  f  (vh+e)wk+1  (4>-<J>h)  (t,xk  +  ||  (t-t  ))dt> 
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For  all  k  €  TL  and  t  €  [tn,  Si+l^ 

|(4>-<Dh)(t,xk)l  <  C  lllflll  (Ax)a(SojJ |+1)  Ax 

so  that 

N  *n+l  i,  n  h 

l  l  !  (vh+e)w£o-<J>  )  (t,xk)dt  < 

n0  Z  tn 

N 

<  C  lllflll  (Ax)a  Z  Z  (l+|ok|)AxAt  < 

no  z 

«£  C  lllflll  ( Ax)  “  I  log  6  I  . 

Consequently 


|q(D+qC2)|  <  c  III  fill  (Ax) a  |  log  6  1 


m  Ilf  lwit+i-wk+1]lv  *EJl*‘*  )ft.*v  ♦  xfCt-t„))dt 


n0  Z  'n 


k  At^  nJ 


Now  from  calculations  analogous  to  the  ones  leading  to  (5.8) 
(5.9),  we  have  in  (and  in  sj!  ) 


|  (v^+e)  ( 4>— 4>^)  |  <  (v^+e)(llij»  II  Ax  +  ll<{>  II  At)  < 

X  °°*Tk  ”»Tk 

<  C(vh+e)  [1  +  ( Vn+  c)  ” 1  +  (vj+e)-1]  Hlflll  (Ax) ' 

k  K 


<  C  lllflll  (Ax) u 


Moreover 


Hence 


n+1  n 
wk  "  wk 


vn+1  vn+1 

1  rvk  "vk  .  vk-l"vku 
8  [  St -  +  - St - ]Ax 


0 "  1  r  n  ,  n  ] 
e  [ok  *  ok.1]  . 


N 

|q(1)+q(2)  |  <  C  lllflll  (Ax) a  |log  6  |  +C  lllflll  (Ax)01  Z  l  |o"|AxAt 


n0  Z 


<  C  lllflll  (Ax)  “  |  log  A  |  . 


This  completes  the  proof  of  Lemma  4.2 


6.  Numerical  Results 


In  this  section  we  discuss  briefly  the  results  of  our 
numerical  experiments  with  the  scheme  (1.11).  All  computa¬ 
tions  were  carried  out  on  the  CDC  6600  at  Indiana  University 
In  each  case  we  specify  only  the  value  of  Ax  used.  The 
values  of  e  and  At  were  always  chosen  to  be  the  smallest 
and  largest  convenient  values,  respectively,  consistent  with 
the  mesh  conditions  [A3]  and  [A4], 

2 

Not  surprisingly,  the  condition  [A4]  on  At/Ax  is 
in  fact  necessary  in  practice.  But  the  condition  [A3]  on  e 
is  probably  overly  restrictive.  For  example,  for  the  speci¬ 
fic  problems  discussed  below,  [A3]  requires  that  e  >  13.8Ax 
which  is  not  "small”  when,  say.  Ax  *  .05  .  Indeed,  for  the 
second  problem  discussed  below,  we  found  that  the  accuracy 
increased  noticeably  as  e  decreased.  A  practical  (but 
not  theoretically  justified)  alternative  to  [A3]  is  the  con¬ 
dition 


(6.1) 


e  *  1 (m-\)  Y0Ax 


Such  a  condition  is  sufficient  for  the  bounds  (2.1)  and 
h  v^ 

(2.2)  for  vn  and  to  remain  in  effect.  (The  more 

stringent  condition  [A3]  was  required  only  for  bound  (2.15) 
92vh 

for  - n-  .)  For  the  specific  problems  discussed  below, 

9x^ 

condition  (6.1)  requires  only  that  e  >  ,58Ax  . 
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For  purposes  of  comparison,  we  used  the  Barenblatt- 
Pattle  solution  v  ,  which  for  m  =  2  is  defined  by 
(see  [16]), 


(6.2) 


TTtyt1  -  c^t)2] 


V(x,t)  =  { 


|x|  <  C ( t) 
I  X I  >  C(t) 


where 


C (t)  =  [12(t+l)]1/3  . 

First,  we  applied  the  scheme  (1.11)  taking  =  v(x^,0) 
»  -C(0)  ,  and  =  c(0)  .  The  computations  were  per¬ 
formed  with  three  different  sets  of  mesh  parameters.  Com¬ 
paring  the  exact  and  computed  solutions,  at  t  =  ^  ,  we 
found  the  following: 


81 


Quite  clearly,  the  observed  errors  in  both  and 

v*1  are  O(Ax)  .  This  is  significantly  better  than  the 

1/12 

rates  predicted  by  Theorem  4.1,  which  are  0(|AxlogAx|  '  ) 

and  0(  |  AxlogAx|  for  c*1  and  v*1  respectively.  These 

discrepancies  are  explained  by  the  fact  that  the  solution 
(6.2")  has  derivatives  of  all  orders  which  are  uniformly 
bounded  on  its  support,  whereas  the  bounds  (4.1)  and  (4.2) 
were  derived  under  the  minimal  smoothness  conditions  which 
all  solutions  are  known  to  satisfy.  Another  difference 
between  the  observed  and  theoretical  results  is  that  the 
observed  rate  of  convergence  for  c*1  is  the  same  as  that 
for  v*1  ,  whereas  Theorem  4.1  predicts  roughly  that 

■  0(1 v-vhH^/2)  .  However,  the  computation  (4.15) 
shows  that,  as  long  as  |c(t)-ch(t)|  is  small, 

U(t)-Ch(t)|  <  r-£—  lv(-,t)-vh(-,t)l  *  C  AX  . 

c2Ct) 

Thus  when  the  interfaces  are  known  to  be  moving  with  speeds 
bounded  away  from  0  ,  the  rate  of  convergence  of  s*1  will 
in  fact  coincide  with  that  of  v  . 

In  the  second  example  we  took  the  same  initial  func¬ 
tion  Vq  as  before,  but  now  with  ^  «*  -3  and  =  3  . 

Thus  Vq  is  neither  concave  nor  continuously  differentiable 
£  r 

on  [Cq  ,  Cq]  .  The  results  were  as  follows: 


We  are  uncertain  as  to  whether  meaningful  comparisons 
can  be  made  between  these  data  and  (4.1)  and  (4.2).  Never¬ 
theless,  it  is  clear  that,  at  least  qualitatively,  Theorem 
4.1  gives  the  correct  result:  in  the  absence  of  smoothness, 
the  convergence  may  be  quite  slow. 

The  scheme  (1.11)  is  thus  seen  to  have  two  shortcomings 
The  first  is  that  the  parabolic  stability  condition  [A4] 
makes  it  impractical  to  apply  the  scheme  with  small  values 
of  Ax  .  This  difficulty  can  probably  be  overcome  by  em¬ 
ploying  instead  a  suitable  implicit  variant  of  (1.11).  We 
intend  to  discuss  such  a  scheme  elsewhere.  The  other  short¬ 
coming  of  the  present  method  is  the  unsatisfactory  rate  of 
convergence.  While  this  phenomenon  is  partly  due  to  the 
weakness  of  the  exact  solutions  themselves,  it  may  be 
possible  to  effect  some  improvement  by  a  more  sophisticated 
treatment  near  the  interfaces. 
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